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Abstract 

We discuss a progress in calculation of Feynman integrals which has been done with help of 
the Gegenbauer Polynomial Technique and demonstrate the results for most complicated parts of 
0(1/A'^"^) contributions to critical exponents of ^''-theory, for any spacetime dimensionality D. 



o 
o 

(N 



> 

(N 
O 



1 Basic Formulae 

Fifteen years ago the method based on the expansion of propagators in Gegenbauer series (see ||l[) 
has been introduced in H, |^. One has shown ||, IJ] that by this method the analytical evaluation of 
counterterms in the minimal subtraction scheme at the 4-loop level in any model and for any composite 
operator was indeed possible. The Gegenbauer Polynomial (GP) technique has been applied successfully 
for propagator-type Feynman diagrams (FD) in many calculations (see ^). In the present paper we 
consider a development of the GP technique (obtained in |^ ) and the application of the results in |^ . 

Throughout the paper we use the following notation. The use of dimensional regularization is as- 
sumed. All the calculations are performed in the space of dimension £> = 4 — 2e. Note that contrary 
to 1^] we analyze FD directly in momentum a;-space which allows us to avoid the appearance of Bessel 
functions. Because we consider here only propagator- type massless FD, we know their dependence on a 
single external momentum beforehand. The point of interest is the coefficient function C/, which depends 
on D = 4 — 2e and is a Laurent series in e. 
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1. First of all, we present useful formulae to use of Gegenbauer polynomials. Following [|| ^, D-space 
integration can be represented in the form 



d'^x 



1 



SD-i{x'Tdx'dx {\ = DI2-1), 



where x — x/^/x^, and Sd-i — 27r'''+^/r(A + 1) is the surface of the unit hypersphere in Ru. The 
Gegenbauer polynomials {t) are defined as j^, ||] 



(l-2rt + r2)-^ = ^ Ciity^ (r<l), C^(l) 



n=0 



T{n + 25) 
n\T{25) '■ 



whence the expansion for the propagator is: 
1 



(a;i - X2) 



25 



X2 



n=0 



{xiT"^' 



Q{xl - x\) + {xl ^ xt) 



(1) 



(2) 



where 



e(y) 



1, ify>0 
0, ify<0 



Orthogonality of the Gegenbauer polynomials C^{x) is expressed by the equation (see [Q) 

A 



C„ [xi £2 ) C„j {X2 X-i ) dx2 



n A 



(3) 



where (5™ is the Kronecker symbol. 
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The following formulae are useful (see [|[ 0): 

~{n - 2py.pW{S) 



p>0 

Substituting the latter equation from for S — X to the first one, we have the following equation 
after the separate analyses at odd and even n: 

^ _ ( ^ ('^ - 2fc + A)r(A) Tjn + S - fc)r(fc + ^ - A) 

- A.^"-2pW ^,p(^) r(n + A + l-A;)r(5-A) ^ ^ 

2. Following ||, |^ we introduce the traceless product (TP) x^'^^■■^" connected with the usual product 
x^'^...x^"- by the following equations 

LLi...ti„ ^ n!(-l)^r(n + A) ^^^^^ „m2p-iM2p ^A'2p+i 

^22Pp!(n-2p)!r(n + A) 



^ (2)2pj5!(n- 2p)!r(n-p + A + 1) ^ ^ ^ 
Comparing Eqs.(^) and (|6|), we obtain the following relations between TP and GP 

2"r(n + A) ^ ^ ^ ' 2"r(2A)r(n + A) ^ ' 



We give also the simple but quite useful conditions: 



(8) 



which follow immediately from the TP definition: ^tii...fj.i...fj.j...fj.„ _ 

The use of the TP x^^'"^^ makes it possible to ignore terms of the type g^^'t^^ that arise upon 
integration: they can be readily recovered from the general structure of the TP. Therefore, in the process 
of integration it is only necessary to follow the coefficient of the leading term x^'^...x'^^. The rules 
to integrate FD containing TP can be found, for example in @, For a loop we have (hereafter 

Dx = {d^x)/{2Tr)^)^. 

1 ..111. ..tin 

jjx - - - - A^-^ia 3) (9) 

j.2a(^^ _ y'^2P (47r)^/2 y2(a+P-\-l) ^ I"' ) 

where 

A (a, (J) = — -— — and a„[a) = — 

an+,n(a + P ~ A - Ij r(Q;) 

Note that in our analysis it is necessary to consider more complicate cases of integration, when the 
integrand contains 6 functions. Indeed, using the Eqs.(||) and (0), we can represent the propagator 
{xi — X2)~'^''^ into the following form Q: 



1 _ ^ 2"r(n + A) „ „ 



L_Q{xl-xl) + (^xj^xtj 



^x,-x,r-^^ nWiX) [xf-'^^^-^^ ~ -^.i (10) 

^The Eg. has been used in |fij, |s|, |lot for calculations of the moments of structure functions of deep inelastic scattering 
^In the case of the propagator (xi — X2)~'^^ with 5 \ we should use also Eq.Ha). 
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Using the GP properties from previous subsection and the connection (Q) between GP and TP, we 
obtain the rules for calculating FD with the 8-terms and TP. 

3. The rules have the following form: 

_ ^ 1 y"" B{m, n|/3, A) 



x2"(a;-y)2/3 ' (4,^)0/2 y2(a+/3-A-l) ^ + + ^ _ A - 1 



(/3=A) 1 V 



(4^)D/2 y2(a-i) r(A) (a- l)(n + A) 

r^2a(^ _ „\2t3 ^^IJ ' (Atv\D/2 „2(a+/3-A-l) 



(/3=A) 1 



^J,l...^J,„ 



(4^)D/2 y2(a-i) r(A) (n + A + l-a)(n + A) 



(11) 



(12) 



where 



B{m,n\P,X) 



r(m + /? + 71) r(m + /3 - A) 



TO!r(TO + n + 1 + A)r(/3) r(/3-A) 

The sum of above diagrams does not contain 0-terms and should reproduce Eq. (|^) . To compare the 
r.h.s. of Eqs.(0[l2|) and the r.h.s. of Eq.(|) we use the transformation of 3-F'2-hypergeometric function 
with unit argument 3i^2(a, b, c;e,b + 1; 1) (see [pl]|): 



E 

k=0 



r(fc + a)r(fc + c) 1 

kW{k + f) k + b 



T{a)T{l - a)r{b)r{c - b) 
T{f-b)T{l + b-a) 

r(i~a)r(a) ^ r(fc + c~/ + i)r(fc + c) i 
r(/-c)r(i + c-/) ^ kir{k + i + c-a) k + c-b ^ ' 



This is the case (when k — m, b = a + (3 — X — 1 , c = n + P) to compare Eq.(^ and the sum of Eqs.(|ll|,|lz 
Analogously to Eqs.(pl]) and (^2|) we have more complicate cases: 



Dx 



1 



2.20(2; „ 2^)2/3 

B{m,n\(3,X) ( (y 



^2{a 



(47r)^/2 



e(y2 - z^) 



y 



.2(q+/3-A-1) 



^2(a+/3-A-i) y\^2j m + a + P-X-1 



TO — a + n+ l + A 



(p=\) 1 



1 



(47r)^/2 r(A) 
1 



1 



6(2/2 - z2) 



?/2(a-i) (a - l)(n + A + 1 - a) 

1 ^e(z2_j/2) ^z2^"+^ e(j/2-z2) 

2 



Da; 



z2("-i)n + AV a-1 
1 



y 



a;2"(a;-?/)2/3 
B(TO,n|/3,A) //y2 



2(a+/3~A-l) 



(/3=A) 1 



e(z^~y^) 
2 " 

/r 

Vz2 



2/' 



All •••Mil 



(47r)^/2 

^ 9(z2-;;2) 

TO + a + Z? — A — 1 



2(a+/3-A-l) 



n + 1 + A — a 
^"'°(«,/3) 



^2^™+/3+» e(y2_^2) 



TO — a + n+ l + A 



(47r)^/2 r(A) 



e(z2-y2) 



2^2(a-l) (Q,_i)(„ + A + l-a) 



1 /e(z2-2/2) /z2x"+^ e(y2-z2) 



z2(a-i)„ + AV a- I 



(2/2) 



rj + 1 + A — Q! 



(14) 



(15) 
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One can easily see that the sum of the above diagrams lead to results identical to 



2 Calculation of complicated FD 

The aim of this section is to demonstrate the result of |^] for a class of master two-loop diagrams containing 
the vertex with two propagators having index 1 or A. 
Consider the following general diagram 



DxDy 



y2"(z-y)2t(^_^)2/32.27(^_y) 



2s 



J{a,t,P,j,s) 



and restrict ourselves to the FD yl(a,/3, 7) — J{a, X,l3,'^, A), which is the one of FD of interest for us 
here. It is easily shown (see [O, |, §) that (cr = 3 + A - (a + /3 + 7)) 



Cf[A{a,l3n)] = Cf[A{a,a,j)] = Cf[J{j,X,X,a,a)] = C f[J{a,j, X, X, 



(16) 



Doing Fourier transformation of both: the diagram A{a, (3, 7) and its solution in the form Cf [A{a, /3, 7)] (z^) 
where hereafter t = X+l — t, t = {a, (3, 7, cr, ...}, and considering the new diagram as one in the momentum 
cc-space we obtain the relation 



Cf[A{a,l3n)] = 



ao(A)ao(a)ao(/3)ao(7) 
ao((5) 



C/[J(5,1,A7,1)] 



= ao(^)"o(a)ao(/3)ao(7)ao(CT) C/[J(/3, l,a,7, 1)] 



(17) 



between the diagram, which contains the vertex with two propagators having the index A, and the similar 
diagram containing the vertex with two propagators having the index 1. 

Repeating the manipulations of e would(see ^, |^) we can obtain the following relations: 



Q[J(A 1, 5, 7, 1)] = Cf[J{fi, 1, a, 7, 1)] - Cf[J{l, 1, 7, d, a)] = Cf[J{d, 1, 1, 7, /?)] 



(18) 



Thus, we have obtained the relations between all diagrams from the class introduced in the beginning of 
this section. Hence, it is necessary to find the solution for one of them. We prefer to analyze the diagram 
A(a,/3,7), that is the content of the next subsection. 

1. We calculate the diagram A{a,(3,^) by the following way^: 



n ^ (10) 2"r(n + A) , 
A(a,/3,7) V )——^ / DxDy- 



n=0 



!r(A) 



Q{z'-y') , Q{y'-z') 



(Qli) 



(47r)^/2 r(A) a 
9(^2 -a;2) 



ri=0 



e(x2 



2"r(n + A) 
n!r(A) 



Dx- 



^tii...fi„^fj.i...fi„ 
a;27(z-a;)2/3 

1 



22(n+A) 
1 



y 



2(n+A) 



A + 71 + 1 — Q! 



(19) 



. 9(z2-x2) 

X + n + a-l' ^22(«+A+a-i) 



e(a;2 - z2) 



x 



2(n+A+Q-l) / 



After some algebra we have got (see |5| ) the result in the form: 

111 



Cf[A{a,Pn)] 



(47r)^ r(A) a-1 



where ^ 



^The symbol = marks the fact that the equation (n) is used on this step. 

''We would like to note that the coefficients in Eqs.(Kh and ( |2l[ ) are similar to ones (see appeared in calculations 

of FD with massive propagators having the mass m. The representation of the results for these diagrams in the form 

!Pn{z^ /m^)" are the coefficients, which are similar to ones in_RqsJ20^ and fclh) is very convenient to obtain the 
results for more complicated FD by integration in respect of m (see f) of results less complicated FD. 



4 



r(n + 2A) 



1 



E 



1 



A^'\a - 1 + 7, /3) + A^'\n + A + 7, /3) 



A + rt + 1 — a 

A^'°{-f, (3) + A"'°{n + a + A + 7-l,/?) 



/ = 



X + n + a — 1 

r(l - f3)T{X + 1 - a)r(A - 1 + a)r(l ~f3 + A)r(l - 7)r(a + /j + 7- A-2) 
r(2A)r(2 + A - a - /3)r(a + 7 - l)r(2 + A - 7 - f3)T{a + /3 ~ A - 1) 
r(n + 2A) (-l)"r(l-/3) 1 



(20) 



E „!r(2A) r(/3-A) \ + n + a~l 
r(a + /3 + 7 - A - 2)r(2 - a - 7) 



(21) 



r(i-7)r(/3 + 7-A-i) 



r(3 - a - /3 - 7 - ?i)r(a + 7 + A - 1 + n) r(7 - A - n)T{2 + 2A - /3 - 7 + ?i) 

Thus, a quite simple solution for A{a,f3,j) is obtained^. In next section we will consider the impor- 
tant special case of these results. 

2. As a simple but important example to apply these results we consider the diagram J(1,1,1,1,q;). 
It arises in the framework of a number of calculations (see ||l2[ |lj, |l^, |2^, |l] ) . Its coefficient function 
I{a) = C/[J(1, 1, 1, 1, a)] can be found (see [||) as follows 



I{a) 



00(1)00(0;) 



C/[J(A,A,A,A,d)] and C/[J(A, A, A, A, d)] = C/[A(5, 3 - A - a, A)] 



ao(a + 2-2A) 
From Eqs. (|2^) and (^ we obtain 

2 r2(A)r(A-a)r(a + l-2A) 



/(a) 



(4^) 



D 



r(2A)r(3A-a - 1) 



(22) 



r2(l/2)r(3A - a - l)r(2A - a)T{a + 1 - 2A) 
r(A)r(2A + 1/2 - a)r(l/2 - 2A + a) 

Note that in |lj] Kazakov has got another result for I {a): 

2 r2(A)r(l-A)r(A-a)r(a + l-2A) 



E 



r(n + 2A) 



^ r(n + a + 1) n + 1 - A + I 



lia) 



(An) 



D 



r(2A)r(a)r(3A-a~l) 



(23) 



r(A)r(2 - A)r(a)r(3A - a - l) 
r(2A- l)r(3- 2A) 



E(- 

ri=0 



yn r(n + 2A) 

' r(n + 2-A) \n + l - X + a 



2A 



From Eqs. (|22| ) and (^3|) we obtain the transformation rule for 3_F'2-hypergeometric function with 
argument —1: 



3^2(20, 6, 1; 6 +1,2- a; -1) = b 
1 - a 



r(2 - a)r{b + a~ l)r{b - a)r(l + a-b) 



r(2a)r(l + - 2a) 



(24) 



b + a-1 



3F2(2a,6,l;6+l,6 + a;l) 



1 + a-b 



3F2{2a, 1 + a - 5, 1; 2 + a - 6, 2 - a; -1), 



where a = X and b = 1 — X + a are used. 

Equation ( p^ ) has been explicitly checked at a = 1 and b = 2 — a (i.e. A = 1 and a = 1), where the 
s^Pj-hypergeometric functions may be calculated exactly. It is very difficult to prove Eq.(p4[) at arbitrary 
a and b values: the general proof seems to be non-trivial. Note that it is different from the equations of 
Ipl] , ^ and may be considered as a new transformation rule. 

^ Before our studies, the possibilitvto represent C'f[A{a,l3,'y)] as a combination of 3F2-hypergeometric functions with 
unit argument, has been observed in 



hv to 
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3 Applications 



The above results have been used for evaluation of very complicated FD which contribute mostly in 
calculations based on various type of 1/N expansions: 

• In the calculation (in |^) of the next-to-leading (NLO) corrections to the value of dynamical mass 
generation (see p^) in the framework of three-dimensional Quantum Electrodynamics. 

• In the evaluation (in ||2^) of the correct value of of the leading order contribution to the /3- function 
of the 0-term in Chern-Simons theory. The /3- function is zero in the framework of usual perturbation 
theory but it takes nonzero values in expansion (see p^). 

• In the evaluation (in ^l[) of NLO corrections to the value of gluon Regge trajectory (see discussions 
in and references therein). 

• In the calculation (in |^) of the next-to-leading corrections to the BFKL intercept of spin- 
dependent part of high-energy asymptotics of hadron-hadron cross-sections. 

• In the calculation (in ||2^) of the next-to- leading corrections to the BFKL equation at arbitrary 
conformal spin. 

• In the evaluation (in |^) of the most complicated parts of 0(1/A^^) contributions to critical expo- 
nents of (/)*-theory, for any spacetime dimensionality D. 

We consider here only basic steps of the last analysis Q . Since the pioneering work of the St Petersburg 
group [ p^ , exploiting conformal invariancc pof of critical phenomena, it was known that the 0{\/N^) 
terms in the large- critical exponents of the non-linear cr-model, or equivalently 0'*-theory, in any number 
D of spacetime dimensions, derives its maximal complexity from a single Feynman integral /(A) (see p^): 



/(A)= Ainn(A,A) 



(25) 



A=0 



where 

y y y2^2(3._y)2A(^_^)2A(y_^)2(A+A) ^ > 

is a two-loop two-point integral, with three dressed propagators, made dimensionless by the appropriate 
power of x^. 



The result, obtained by GP technique, is 



/(A) = .1.(1) - .1.(1 - A) + > (27) 



where ^{x) = T'{x)/r{x) and 

/•I ^2A-l 

^{X)^A dx- j{Li2(-a;)-Li2(-l)} (Li2(x) = V xVn^), (28) 

In 1^ the integral /(A) has been expanded near D — 2 and D = 3, respectively, (i.e. for D = 2 — 2e 
and D = 3 — 2e) up to in the form of alternative and non-alternative double Euler sums 
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